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ABSTRACT 

o 

Context. The existence of dusty debris disks around a large fraction of solar type main-sequence stars, inferred from excess far-IR 
and submillimetre emission compared to that expected from stellar photospheres, suggests that leftover planetesimal belts analogous 
Cw to the asteroid- and comet reservoirs of the Solar System are common. 

Aims. In order to detect and characterise cold extended dust originating from collisions of small bodies in disks, belts, or rings at 
Kuiper-Belt distances (30-50 AU or beyond) sensitive submillimetre observations are essential. Measurements of the flux densities 
at these wavelengths will extend existing IR photometry and permit more detailed modelling of the Rayleigh- Jeans tail of the disks 
spectral energy distribution (SED), effectively constraining dust properties and disk extensions. By observing stars spanning from a 
few up to several hundred Myr, the evolution of debris disks during crucial phases of planet formation can be studied. 
Methods. We have performed 870yimr observations of 22 exo-Kuiper-Belt candidates, as part of a Large Programme with the 
LABOCA bolometer at the APEX telescope. Dust masses (or upper limits) were calculated from integrated 870yum fluxes, and 
fits to the SED of detected sources revealed the fractional dust luminosities / dust , dust temperatures r dllst , and power-law exponents f3 

l/~) of the opacity law. 

Results. A total of 10 detections with at least 3tr significance were made, out of which five (HD 95086, HD 131835, HD 161868, 
HD 170773, and HD 207129) have previously never been detected at submillimetre wavelengths. Three additional sources are 

t-H marginally detected with > 2.5cr significance. The best-fit ji parameters all lie between 0.1 and 0.8, in agreement with previous 

results indicating the presence of grains that are significantly larger than those in the ISM. From our relatively small sample we es- 
timate / dllst x r" , with a ~ 0.8-2.0, and identify an evolution of the characteristic radial dust distance R dust that is consistent with the 
t I/3 increase predicted from models of self-stirred collisions in debris disks. 



Key words. Stars: circumstellar matter - Stars: individual: HD 105, HD 21997, HD 30447, HD 95086, HD 98800, HD 109573, 
HD 131835, HD 141569, HD 152404, HD 161868, HD 170773, HD 195627, HD 207129 - Stars: planetary systems: formation - Stars: 
planetary systems: planetary disks 

1. Introduction stellar photospheric contribution) discovered by the Infrared 



Astronomical Satellite (IRAS, e.g. Aumann et al. 



The existence of warm dust around solar type main-sequence [Neugebauer etlLl[T984l, and studied by the Infrared Space 



1984 



stars is implied from infrared (IR) excess emission (above the 0bservatory (ISO, e.g. |Kessler et al.|[T996|, and Spitzer (e.g 



[Werner et al.||2004] l. Since such dust has a limited lifetime 
Send offprint requests to: R. Nilsson its existence suggests the presence of larger asteroidal and/or 

'Based on observations with APEX, Llano Chajnantor, Chile (OSO cometary bodies which through collisions continuously replen- 
programme 081.F-9330(A)). 
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ish the dusty debris disk. Excess emission at wavelengths of 
25 /mi and shorter comes predominantly from hot dust (>300 K), 
located within 1-5 AU of the star (similar to the asteroidal belt in 
our Solar system), while an excess at 25-100 /mi probes cooler 
dust (~ 100-300 K) located at 5-10 AU. However, in order to 
detect cold and very extended dust, originating in Kuiper-Belt 
analogs (on the order of 100 AU from the star), observations at 
submillimetre (submm) wavelengths are needed. Although very 
few older stars show evidence of hot dust, likely due to clear- 
ing by planets (e.g. Beichma n et al.|20 06b, Me yer et al.|2 008), 
cool debris disks appear more common and many should be de- 
tectable at submm wavelengths based on extrapolation from ob- 
served far-IR fluxes. 

The potential correlation of cold dust and planetesimal belts 
with the formation of exo-planets in these systems is another 
incentive for larger submm surveys of nearby solar-type main- 
sequence stars. Beichman et al. ( 2006a I and Bryden et al. ( 2006 ) 
reported Spitzer MIPS observations of a large sample of nearby 
F, G, and K main-sequence stars, and detected an overall ex- 
cess of 70/mi emission towards 13% of the stars. From a larger 
sample of solar-type stars, spanning ages between 3 Myr and 
3 Gyr, studied in the Spitzer Legacy programme Formation and 
Evolution of Planetary Systems (FEPS), showed a total 70/im 
excess rate of 7% with high (50 times the photosphere) frac- 
tional excesses at ages between 30 and 200 Myr, and signifi- 
cantly lower fractional excess for older stars (Carpenter et al. 
|2009l|Hlfle"nbrand et al.|20081 >. There was no clear tre nd for the 
temporal evolution of 70 /mi excess in FEPS, however Su et al. 
(2006) estimate a decay time of ~400Myr for A type stars. At 
submm and mm wavelengths 3% of the FEPS stars show evi- 
dence of continuum emission, with a suggested decrease in dust 
masses and/or changes in the grain properties for stellar ages of 
10-30 Myr ( |Carpenter et alipOOB) . 

An additional advantage in making observations in the 
submm region is that measured integrated fluxes are directly pro- 
portional to the temperature and mass of the disk, due to the fact 
that they (in most cases) sample the Rayleigh-Jeans tail of the 
spectral energy distribution (SED), and that the disk can be as- 
sumed to be optically thin at these wavelengths. At the same 
time, the thermal radiation in the submm is dominated by large 
and cool dust grains, thus providing a better estimate of the total 
dust mass (although still giving just a lower limit). By combining 
submm observations with IR-photometry a fit to the SED can be 
made, yielding constraints on the temperature and radial extent 
of spatially unresolved disks. 

In this paper we present the first results of a Large 
Programme submm-survey with the Large APEX Bolometer 
CAmera (LABOCA) at the Atacama Pathfinder Experiment 
(APEX) telescope, targeting stars of spectral type B to M, be- 
longing to nearby stellar associations with ages ranging between 
10 and 100 Myr. While the main aim of the study is to investigate 
disk evolution by relating disk parameters and grain properties 
to fundamental stellar parameters (metallicity, mass, and age), 
we also hope to spatially resolve extended dust disks around 
some of the most nearby stars. The observational scheme is de- 
signed to reach ~2 mJy/beam root-mean-square (RMS) sensitiv- 
ity, enabling the detection of debris disk masses as low as the 
mass of the Moon (comparable to the Solar System's Kuiper- 
Belt mass of 3 MMoon; Iorio 2007 and references therein), and to 
make observations of exo-Kuiper-Belts out to distances of about 
50-100pc. A successful precursor study of the ~12Myr-old /3 
Pictoris Moving Group clearly detected 2 out of 7 stars (Nilsson 
|et al.|2009) . 



2. Stellar Sample 

The sample of 22 F-, G-, and K-type stars selected for this 
first round of observations were based on levels of far-IR 
excess (measured by IRAS, ISO and Spitzer) and proximity 
(also making them targets for Herschel Space Observatory Key 
Programmes on debris disks), in order to optimise initial detec- 
tion rates and get meaningful upper limits. Positions, spectral 
types, ages, and distances of the objects are listed in TablefT] 
together with integration times of the observations. The age es- 
timates of most stars, which are based on lithium abundances, 
H-R diagram location, moving group (MG) association, chromo- 
spheric activity, etc., all carry large uncertainties. E.g., HD 105 
is a member of the Tuc-Hor MG (Mamajek et al. 2004), with 
an estimated age of 27+ 1 1 Myr ( |Mentuch et al.|20 08 ), but could 
be as old as 225 Myr based on the level of chromospheric activ- 
ity ( |Apai et al.|[2058l ; and HD 98800, HD 25457, HD 141569, 
HD 152404, and HD 207 129 might still be in their pre-main- 
sequence phase. Due to these inconsistencies we have adopted 
probable approximate ages with upper and lower limits that en- 
compass the range of published values from different methods. 
Consequently, our ability to resolve details about the temporal 
evolution of disk properties in the young 10-100 Myr systems, 
where terrestrial planet formation could still be in progress, ac- 
cording to, e.g., isotopic dating of meteorites and the Earth- 
Moon system (Touboul et al. 2007 1 and dynamical modelling 
( |Mandell et al.|2007) , is somewhat restricted. 



3. Observations and Data Reduction 



The 12-m diameter submm telescope APEX (Glisten et al. 
2006) is located at an altitude of 5100m in the Chilean Andes, 
offering excellent atmospheric transparency. It is operated 
jointly by Onsala Space Observatory, the Max-Planck-Institut 
fur Radioastronomie, and the European Southern Observatory. 
Since 2007, the LABOCA bolometer array has been available for 
observations at a central wavelength of 870/mi (345 GHz) and 
bandwidth of 150/mi (60GHz), covering a 1L4 field-of-view 
with its array of 295 bolometers; each with a nearly circular 1972 



full-width at half-maximum (FWHM) Gaussian beam ( Siringo 
|et al.|[2009] ). A spiral pattern mapping mode was employed in 
order to recover fully sampled maps from the under-sampled 
bolometer array. Siringo et al. (2009) provide a detailed descrip- 
tion of the instrument, observing modes, performance and sen- 
sitivity. 

Observations of the targets presented in Table [T] were ob- 
tained between June 4 and October 22, 2008, amounting to 
roughly 80 out of 200 allocated hours of this Large Programme. 
Several individual 7.5 min long spiral scans, producing raw maps 
with an approximately uniform noise distribution within a ra- 
dius of about 4' around the central position of the source, were 
obtained for each star. In-between these scans, skydips (to de- 
termine the correction for the atmospheric opacity), flux calibra- 
tion, and measurements of pointing accuracy and focussing were 
performed on selected calibration objects. 

The scans were reduced and combined using MiniCRUSH 
(vl.05), an adapted version of CRUSH (Comprehensive 
Reduction Utility for SHARC-2, Kovacs 2008), with further pro- 
cessing for analysis and plotting of final maps in Matlab (v7.9.0). 
As all of the sources were relatively faint we used filtering and 
settings optimised for point source extraction, with full beam 
smoothing, resulting in clean maps with 27" effective resolu- 
tion. Details of the data reduction procedure can be found in 
Kovacs (2008) and are summarised in Nilsson et al. (2009). The 
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Table 1. Observing log with stellar properties and integration times. 



Object ID Name 



R.A." 
(h m s) 



Dec. 

(° ' ") 



Spectral Type 



— T upper 

» lower 

(Myr) 



Reference 



30^ 
300™ 

-"'21)0 

50 150 
2Q ioo 

13004000 

18of° 

220 320 
20 50 

wf 

10° 
14I8 

1H 10 

200 300 
5 8 

n)l 

180 



Distance"" 
(pc) 



Integration time 
(min) 



HD105 
HD 17390 
HD 21997 
HD 25457 
HD 30447 
HD 31392 
HD 32297 
HD 61005 
HD 78702 
HD 95086 
HD 98800 
HD 109573 
HD 131835 
HD 139664 
HD 141569 
HD 152404 
HD 161868 
HD 164249 
HD 170773 
HD 193307 
HD 195627 
HD 207 129 



TVCrt 
HR4796 



AKSco 
yOph 



00 05 

02 46 

03 31 

04 02 
04 46 

04 54 

05 02 
07 35 

09 09 

10 57 
1122 
12 36 

14 56 

15 41 

15 49 

16 54 

17 47 

18 03 
18 33 
20 21 
20 35 
2148 



52.54 
45.10 
53.65 
36.74 
49.53 
04.21 
27.44 
47.46 
04.21 
03.02 
05.29 
01.03 
54.47 
11.38 
57.75 
44.85 
53.56 
03.41 
00.92 
41.03 
34.85 
15.75 



-5145 11.0 
-21 38 22.3 
-25 36 50.9 
-00 16 08.1 
-26 18 08.8 
-35 24 16.3 
+07 27 39.7 
-32 12 14.0 
-18 19 42.8 
-68 40 02.4 
-24 46 39.8 
-39 52 10.22 
-35 4143.6 
-44 39 40.3 
-03 55 16.4 
-36 53 18.6 
+02 42 26.2 
-5138 56.4 
-39 53 31.3 
-49 59 57.9 
-60 34 54.3 
-47 18 13.0 



GOV 
F3IV/V 
A3IV/V 

F5V 

F3V 

K0V 

A0V 

G8V 

A0/A1V 

A8III 

K5V+K7V+M1V+?" 

A0V+M2.5 / 

A2IV 

F4V 

A0Ve+M2+M4 s 

F5V+F5V' 1 

A0V 

F5V 

F5V 

GOV 

F0V 

G2V 



id 

170 

12? 

501) 



200 
5700*950 

5100 
30f° 

3000fij 00 



1,2,27 

3,4 

5,3,6 

7 

6,8 

9 

10 

11 

6 

12, 13, 14 

15 

16, 17, 18 

8,6 

7 

19,26 

20 

16, 17 

21 

3,5 

22,27 

3, 12 

. 22, 24, 25, 27 



40 
45 
74 
19 
78 
26 
112 
35 
80 
92 
47 
67 
111 
18 
99 
145 
29 
47 
36 
32 
28 
16 



104 
24 
52 
97 
65 
53 
39 
53 
91 
32 
62 
185 
145 
170 
206 
127 
218 
166 
311 
90 
98 
354 



References. (1) Hollenbach et al. (2005) and references therein; (2) Apai et al. (2008); (3) Zuckerman & Song (2004b); (4) Zuckerman & Song 
( |2004a| ); (5 )|Song et al.H2003 1; (6) Mo oret al.| ( |2006| >; (7) |L6pez-Santiago et al.|l|2006) and references therein; (8)|Rhee et al. ( 200 7[> and reference s 
ther ein; (9) Carpenter et al. ( 2009 ) an d references therein ; (lO) IPerryman et al.|(|1997T; (1 l)|Roccatagliata et al.|<|2009) >; (12 )|Lowrance et al.|<|2000| ; 
(13)|deZeeuwet al. (1999); (14) Mamajek et al. (2002); (15)|Barrado Y Navascue s (2006); (16) Song et al. (20 01]; (17)[Rieke et al. (2005 ); (18) 
Stauffer et al.|(|1995| ); (19)|Weinberger et al. (2000); ( 20)|Chenet al.|(|20 05) and references therein; (21) Zuckerman et al. (2001 ); (22) Mamajek & 
Hillenb rand (20 08J; (23)|Henry et al.| l |1996V ; (24) |Bryden et al.| ( |2006| > and references therein; (25) |Zuckerman & WebblfeoOO) ; (26) |Merfn eTaT 
( 12004) ; (27) |Valenti & FischerH2005) . 



Notes. <°> Positions are from the Hipparcos Catalogue ( [Perryman et al.|l997 1. <6) Spectral types are from Torres et al. (2006); Thoren & Feltzing 
(2000); Manoj et al. (2006); Th evenin & I diart ( 1999}. ,CJ Distances are from the Hipparcos Catalogue (Perryman et al.|1997| > w From Lopez- 
Santiago et al. (2006). Marked as T Tau-type in SIMBAD astronomical database. w Quadruple system composed of two spectroscopic binaries 
Aa+b and Ba+b, separated by 0/848 iTPrato et al. ||2001 1, IR excess originates from the Ba+b (K7+M1V binary) component (for which the 
coordinates are given) (see Laskar et al.||2009| and references therein). tf) Binary system with the secondary (M2.5 pre-main-sequence) star 
located 7/7 to the s outhwest of t he primary (A0V). IR excess comes from HD 109573A (also designated HR4796A or GQLup) with an imaged 
circumstellar ring l Schneider et al.|[200 9). (s) Primary component HD 141569 A, which is possibly still in its pre-main-sequence stage, has a 
binary companion consisting of two M stars at ~ 8" separation. The double pair could be unbound and merely flying by (Reche et al. 2009). 
"" Spectroscopic binary with high IR excess attributed to long-lived primordial, circumbinary disk (Uzpen et al. 2007 , Gomez de Castro 2009). 



general reduction steps involve flux calibration (opacity correc- 
tion and counts-to-Jy conversion), flagging of bad (unresponsive, 
very noisy, or those that were too fast or too slow) channels, cor- 
related noise removal, despiking, data weighting, and map mak- 
ing. Siringo et al. (2009) can be consulted for a description of 
these steps and their relation to the observing mode of the in- 
strument. 

The final maps had a sensitivity between 2 and 10 mJy/beam. 
The integrated flux density of detected sources was found by fit- 
ting a 2-D Gaussian to the source region after baseline subtrac- 
tion and gradient correction, with errors estimated from root- 
mean-square (RMS) noise calculation in the specific integration 
region, together with an absolute calibration error estimated to 
be 10% ( [Siringo et aLpOO^ . 



4. Results 

Of the 22 far-IR excess stars observed, 10 (HD 21997, 
HD95086, HD98800, HD 109573, HD 131835, HD 141569, 
HD 152404, HD 161868, HD 170773, and HD 207129) were de- 
tected with at least a 3cr significance, giving a detection rate 



of -45%. Five of the debris disks (HD 95086, HD 131835, 
HD 161868, HD 170773, and HD 207129) have previously never 
been seen at submm wavelengths. HD 105, HD 30447, and 
HD 195627 are listed as marginally detected in Table|2] show- 
ing flux peaks of > 2.5cr significance at the position of the star. 

The final 870 /mi flux-density maps of the detected and 
marginally detected objects are found in Fig.fT] which show the 
flux level color-coded in a scale ranging from zero to the peak 
flux density (in Jy/beam) of each individual map. The dotted 
contour represents the lcr-level and the first solid contour out- 
lines the 2cr-level, with subsequent contour lines spaced by lcr. 
A circle representing the effective 27" beam size after smooth- 
ing, and a line showing the angular size subtended by 1000 AU 
at the distance of the object, has been inserted in the lower left 
comer. 

Several significant flux density peaks can be found sur- 
rounding the position of the observed star, also in parts of the 
LABOCA field-of-view outside the plotted 200 x 200" maps. 
Considering the (down to) 2mJy/beam sensitivity we are ex- 
pected to hit the extragalactic confusion limit, permitting de- 
tection of many background submm galaxies (see e.g., Wang 
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Table 2. Integrated flux density, root-mean-square noise levels, and derived dust temperature, power-law exponent of the opacity 
law, mass, fractional dust luminosity, and characteristic radial dust distances for the 22 stars observed at 870-yum. 



Object ID 


Integrated Flux Density, F 


RMS Noise, a 


r r a 

1 dust 


P 


Dust Mass'' , M dust 


./dust 


^dust 




(mJy) 


(mly/beam) 


(K) 




(™Moon) 


(10-*) 


(AU) 








DETECTIONS (> 3cr) 










HD 21997 


17.6±8.0 


4.6 




52 


0.7 


33±L5~ 


4.9 


240 


HD 95086 


41.3+18.4 


10.4 




150+(15-45) t ' 




>76 J 


6.4-14 


16+(250-1500) 


HD 98800 


33.6±8.4 


3.9 




155" 


0.1 


8.5 


470 


2.4 


HD 109573 


21.5+6.6 


3.1 




88^' 


0.7 


19.5±6.0 


38 


77 


HD 131835 


8.5±4.4 


3.0 




89 


0.1 


21+11 


25 


47 


HD 141569 


12.6±4.6 


2.3 




80« 


0.8 


27±10 


73 


110 


HD 152404 


42.9±9.8 


3.7 




— s — 










HD 161868 


12.8+5.2 


2.5 




73 


0.5 


2.6+1.1 


0.90 


120 


HD 1707/3 


18.0±5.4 


2.6 




43 


0.7 


9.7±2.9 


4.5 


170 


HD 207129 


5.1±2.7 


2.1 




38 


0.7 


0.6±0.3 


0.94 


130 


MARGINAL DETECTIONS (> 2.5<r) 


HD105 


10.7±5.9 


3.7 




40 


0.6 


7.7±4.2 


2.6 


110 


HD 30447 


6.9±5.0 


4.5 




61 


0.3 


12.3±8.9 


9.3 


56 


HD 195627 


13.0+.7.I 


4.1 




57 


0.3 


3.1+1.7 


1.0 


74 


NON-DETECTIONS 


HD 17390 




6.8 




55' 




<13.4 


W 




HD 25457 




3.3 




70' 




<0.9 


1.0 s 




HD 31392 




4.7 




49' 




<3.4 


1.6* 




HD 32297 




6.5 




30' 




<150 


33^ 




HD 61005 




6.0 




58* 




<6.6 


25* 




HD 78702 




4.2 




36'" 


2.0"' 


<39 


2.6'" 




HD 139664 




3.5 




78" 




<0.7 


1.2" 




HD 164249 




3.5 




78 




<5.3 


5.9 U 




HD I933UV 




3.8 















Notes. I^JFor detected and marginally detected sources we use the dust temperature derived from the best fit to the spectral energy distribution. 
EMlUpper 3<r limit on the dust mass for undetected sources. <c) In order not to be resolved the second dust belt should be at <1500AU, which 
would mean >15K, while a temperature > 45K would not be compatible with far-IR data. W) Calculated for upper limit on dust temperature. 
EH Ambiguou s fit (see discussion in S ectionlS^. iH Source appears spatially res olved (see discussion in Appendix [Ap. See also more detailed 
modellin g by Augereau et al. 1 1999b| and SCUBA me asurement by Sheret et al. 1 2004 i.lSJ Source appears spat ially resolved ( see discussion in 
Appendix|A|(. See also modelling by Sheret et al. (2004). " ;) No simple debris disk model can be fitted (see Section p^ , (,) From Rhee et al. (2007). 
El Fro m Moor et al. (2006). ( *> From Hillenbra rdeTSl ( |2008| ).ESFrom [Maness et al.| [2008} . E] From |Williams & Andrews| l l2006fr .lMFrorn 
Beichman etal.| {2006a). (o) From Re bull et al"| |2"008). <p) Not enough photometry data for reliable SED fit. 



|et al. 2004| ). This also effects the interpretation of our results 
for HD 109573 and HD 141569, which appear spatially resolved 
even with our rather large effective beam which had been opti- 
mised for point source extraction in the reductions. These issues 
are discussed further in Appendix [A] 

The measured integrated fluxes and RMS noise levels are 
listed in TableE] together with the disk parameters (temperature, 
power-law exponent of the opacity law, dust mass, fractional 
dust luminosity, and characteristic radial distance of the dust) 
derived from the analysis outlined in Section|5] For undetected 
sources, an upper 3<x limit on the dust mass is given. 



5. Discussion 



5.1. Frequency of Cold Dust Disks 



Similar to the results of Carpent er et al.| ( 2005| ), we find that 
about half of the stars previously observed to have far-IR fluxes 
in excess of stellar photospheric flux also have an excess of 
submm continuum emission. The five new submm disks dis- 
covered increase the number of known exo-Kuiper-Belt systems 
from 29 to 34 (see Appendix fBl, counting the possible tran- 
sitional disk£] (HD 98800, HD 141569, HD 152404), and with 



'HD 98800, HD 14 1569, and HD 15240 4 are, according to sp ectro- 
scopic studies by e.g. |Furlan et al.| (2007), |Merfn et al.| (|2004[), and 



three additional sources (HD 105, HD 30447, and HD 195627) 
awaiting confirmation. Although it is hard to estimate the total 
frequency of cold extended disks around main-sequence stars, 
it is likely to be at least the estimate of 10-15% suggested 
from IRAS ( e.g. |Backman & Parescel[T993| |Rhee et al.||2007] l 
and Spitzer (Hillenbrand et al. 2008 ) detections of debris disks 
around main-sequence stars, and may be up to 25% including 
disks too faint or cold to have been detected with current instru- 
ment sensitivities (Matthews et al.|2007) . 

An interesting result is that all four multiple systems in our 
sample, HD 98800, HD 109573, HD 141569, and HD 152404, 
were detected at 870 fim. Multiple stars have previously been 
deliberately left out from most debris disk surveys under the as- 
sumption that the formation of planetesimals in such systems 
would be inhibited. However, the Spitzer MIPS observations pre- 
sented by Trilling et al. (2007) show that the incidence of debris 
disks around main-sequence A3-F8 wide (>50AU) binaries is 
comparable to that around single stars and is even marginally 
higher for tight (<3 AU separation) binaries. Our results con- 
firm these trends, since the four multiple star systems in our 
sample are three wide binaries and one spectroscopic tight bi- 



Alen car et al.| 12003 1, respectively, pre-main-sequence stars and prob- 
ably do not possess actual debris disks in the sense that they contains 
a large amount of gas, and might thus fall within the loosely defined 
transition disk category. 
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nary (HD 152404). These results should however be taken with 
caution because of the morphological states of these four mul- 
tiple stars systems. The main one is that all four systems are 
young (~5— 15 Myr old) and thus should naturally have a higher 
probability of possessing massive debris disks (see discussion 
in Section 5.5 1. This bias towards young stars in our sample of 



multiple star systems make it difficult to infer trends related to 
the binarity. 

The fraction of detected debris disks in this first round of 
observations (45%) is considerably higher than that of previ- 
ous submm surveys, e.g. that of nearby bright stars (9%, Holmes 
etal.|2003),)SPictoris Moving Group (17%,[Nilsson et al.|2009 



Liu et al.|2004) , FEPS nearby stars (3%, |Carpenter et al.|2005) , 
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nearby G stars (15% , |Greaves et al.|2005| l, and M dwarfs (6%, 
|Lestrade et al.|2006) . This cannot only be explained by higher 
instrument sensitivity but is most likely due to a selection ef- 
fect, having chosen an initial sample of stars with high far-IR 
excesses, and we are not expecting such a high detection rate 
from our ongoing Large Programme observations. 



As anticipated with RMS noise levels as low as 2 mJy/beam 
we find a number of significant flux density peaks originat- 
ing from background submm galaxies. Based on results from 
previous deep galaxy surveys at mm and submm wavelengths 
(e.g. |WeiB et al . 2009 1 Bertoldi et aLp007l |Ivison et al-pOUT] 



Clements et al.|2004[l we would expect about 4-6 extragalactic 



sources to be detected (at >3cr levels) in the LABOCA field- 
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Fig. 1. Maps of objects detected at 870//m with LABOCA at the APEX telescope. Each map has been smoothed with a circular 
Gaussian corresponding to the FWHM of the beam, and the resulting effective FWHM image resolution of 27" is shown as a circle 
in the maps. In addition, a scale bar corresponding to 1000 AU at the distance of the star is inserted. The first solid contour represents 
2cr flux levels, with the following contours at increments of lcr (or 2<x for HD 98800 and HD 152404), while the dotted contour 
outlines the lcr level. The position of the star is at (Aa, Ac?) = (0", 0"). 



of-view with this sensitivity. This is consistent with the average 
number of four 3<x peaks within a 4' radius around the central 
position of the source that we observe in our final maps. We 
can estimate the probability of a finding a random background 
source within an angular distance r of a given position f rom the 
expression P(< r) = 1 — exp(npr 2 ) ( jCondon et al.|l998 1, where 
p is the source density per square degree. With a 27" effective 
beam we calculate a roughly 2% likelihood of a chance align- 
ment. Although the risk of a background galaxy falling within 
the central beam is so small, the probability that it does must be 
considered, as it could produce a false disk detection or a seem- 
ingly extended disk. With the commissioning of the Atacama 
Large Millimetre/Submillimetre Array (ALMA) now on its way, 
unprecedented angular resolution (>0701 , Peck & Beasley 2008 ) 
at submm and mm wavelengths (from 350 /zm to 10 mm) will 
soon permit such ambiguities to be resolved. 

5.2. Spectral Energy Distributions 

By combining the integrated submm flux of our detected objects 
with photometry from optical and infrared measurements found 
in the literature, we made fits to the SED in order to find the dust 
temperature and radial extent of the disk. The star is modelled as 
a simple scaled blackbody, while the disk is approximated with 
a modified scaled blackbody function, as 



F disk(v) = 



2hv 3 






kqD. 



(1) 



where h is Planck's constant, k is the Boltzmann constant, and 
c is the speed of light, with the radiating surface area and dis- 
tance incorporated in the scale factor C. Here we assume that 
the opacity index varies with the frequency as a power law, i.e. 
k v = K^ivlvof, and that the grains dominating the flux in the IR 
and submm all have the same temperature. The latter simplifying 
assumption is applied since we do not have a spatially resolved 



source at submm wavelengths, although it would be reasonable 
to believe that colder dust located at Kuiper-belt distances would 
be responsible for most of the submm emission, while warmer 
dust located at asteroid-belt distances would dominate in the far- 
IR. Best-fit parameters were found by minimizing the reduced 
X 2 , i.e. the sum of squared errors weighted with the accuracy of 
the individual photometry measurements, divided by the num- 
ber of degrees of freedom. Resulting SED fits are presented in 
Fig. [2] while values of the dust temperature, rd ust , and the power 
law exponent of the opacity law, /?, for the 10 detected (and 3 
marginally detected) objects are listed in Table|2] 

Although only the best-fit /3 parameter is quoted (without er- 
rors), we note that in almost all fits a value between 0.2 and 
0.7 gives roughly the same x 1 - The difficulty in determining f3 
is due to relatively large errors in measured submm (and longer 
wavelength) data, in addition to the possible existence of colder 
dust that would have been better modelled with a second dust 
component, but now instead shifts /3 in a single component fit 
towards zero. One example is HD 95086 which clearly cannot 
be modelled with a single component dust disk. For this star 
we added a second disk with a lower temperature limited by 
the maximum extent of the unresolved disk (<1500AU) and 
upper limit constrained by the two IR points at 13 and 33 //m 
( Chen et al.|2006) . The fit shown in Fig. |2(e) is for a lower limit 
5 K disk. Overall, the derived /3 values for our sources indicate 
significantly larger (and perhaps amorphous and/or fractal) dust 
grains (Miyake & Nakagawa 1993 1 Mannings & Emerson 1994 



Pollack et al. 1994 1 in circumstellar debris disks compared to 



the mostly unprocessed grains in the interstellar medium (ISM), 
which have /3~2 ( |Hildebrand|[T983| |Li & Draine"|20"0"T]>. This is 
in agreement with previous results by e.g. Najita & Williams 
( |2003] l, [Niis"son et al.|p009] l, and |Roccatagliata et al.| ( |2009| l who 
found j8=0-l for various debris disks. 

In three cases the SED modelling does not yield a good fit. 



For the quadruple system HD 98800 (Fig. 2(e) i previous submm 



and mm photometry has produced somewhat differing results, 
viz. the 1350jt/m flux from Sylvester et al. (2001 ) is significantly 
higher than the 1300/im flux from Sylvester et al.| ([T996), with 
our 870 /im measurement more consistent with the latter for a 
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typical Rayleigh-Jeans slope. It can not be excluded that pro- 
cesses related to the orbital phase of the tight (1 AU) binary 
HD 98800B cause this seemingly bimodal flux variation of its 
circumbinary disk. Photometric variability has been found at 
shorter wavelengths (Soderblom et al. 1998}. In any case, this 



unusual system of two eccentric binaries orbiting each other with 
highly inclined orbital planes (see e.g. Verrier & Evans 2008 and 



references therein) should have a complicated dynamical effect 
on the debris (or transitional) disk around the B pair, e.g. the high 
IR-excess has been proposed to originate from a puffed up outer 
dust ring which is gravitationally perturbed by the HD98800A 
binary at perihelion (50 AU) passage (Furlan et al. 2007). 
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For HD 141569 the mid-IR flux is difficult to fit consistently 
with longer wavelength data in a simple disk model, however 
our results are very similar to those of Sheret et al. (2004). 

The photometry data of HD 152404 (AKSco) can clearly 
not be fitted with a simple modified blackbody emitting disk. 
In Fig. |2(i)| we have inserted a dashed line representing the SED 
of the binary consisting of two F5V stars, showing that excess 



emission is present already at A<1 /im. This is indicative of hot 
dust in the system's young circumbinary disk, and it is probably 
still in its protoplanetary or transitional disk stage. 

Although more detailed modelling of the SED could have 
been performed (and for some sources previously has been made 
by other authors) the lack of sensitive photometry at wave- 
lengths between ~1 00-500 /mi still hampers our ability to accu- 



R. Nilsson et al.: Kuiper Belts Around Nearby Stars 




100 1000 

Wavelength (^im) 



(111) 

Fig. 2. The total modelled SEDs of detected objects are shown as a solid black line, with the stellar and disk components represented 
by red dashed and dotted lines, respectively. Our obtained 870 /im integrated flux density is given by a green square, while blue 
circular markers (and triangles symbolizing upper limits) are data from the literature. Fits which yield exceptionally large x 1 values 
and other special cases are discussed in Section |5T2"1 The 800 pm flux measured for HD 141569 (red marker) at CSO (|Walker & 
Butner 1995) has been left out of the fit since it is probably incorrect. References: HD 105 (Moor et al. 2006; Hillenbrand et al. 
2008 ; Carpenter et al. 2005 ; Cutri etal.|2003"l); HD 21997 (Moor et al. 2006; Williams & Andrews|2006 ; Cutri et al. 2003 ); HD 30447 



(Moor et al. 2006; Cutri et al. 2003); HD 95086 (Chen et al. 2006; Cutri et al. 2003); HD 98800 (Low et al. 2005 ; Schiitz et al. 2005 
Mannings & Barlow|1998[ |Cutri et al.|2003||Walker & Heinricnsen|2000[ |Gregorio-Hetem et al.|1992||Garcia-Lario et al.p990 



Sylvester et al.|1996f|Cieza|2008|ISylvester et al.|2001|>; HP 109573 (|Mo6r et al.|2006[|Sheret et al.|26()4||Barrado Y Navascues 



2006 , [Cutri et al.|2003[|Low et al.|2005||Jura et al.|19 95; Teles co et al.|2000[|Jura et"aT. 



'993 T JFajai-do-A costa et al.|199 8; Wahhaj 
2001 1996; Beichman et al. 1988 ; Walker 



et al.|2005|l; HD 131835 (Mooret al. 2006| |Cutri et al.|2003|l; HD 141569 (Sylveste r etal. 



& Butner 1995 ; Sheret et al. 2004); HD 152404 (Chen et al. 2005 Jensen et al. 1996 ; Acke et al. 2004; Cutri et al. 2003 ); HD 161868 
([Suet al. 2008; Fajardo-Aco sta et al.|1999"l|Jaschek et al.|1991||Cutri et al.|2003); HP 170773 (JMoor et al 12 006; C utri et al.|2003|l; 



HD 195627 (|Rebull et al.|200 8 ; Sylvester & Mannings 2000, Mannings & Barlow|199"8"l|Cutri et al.|2003[); HP 20712 9 ([Mannings 
& Barlow|1998[|Trilling et al.|2008| |Schiitz et al.|2005| |Cutri et al.|2003||Mamajek et al.|2004[|Laureijs et al.||2002| |Jourdain de 
Muizon et al.|1999[|Tanner et al.|2009| )~ 



rately determine grain properties and size distributions in cold 
debris disks. This will be addressed with ongoing Herschel 
PACS/SPIRE observations of disks, which are filling the gap 
between current far-IR and submm measurements, and will de- 
tect dust masses down to a few Kuiper-Belt masses for nearby 
stars ( Panchi et al.||2010 l. With such data, any deviation from 
the theoretical power law size distribution of observable dust, 
e.g. the "wavy" size distribution predicted from numerical mod- 
els of collisional processes in debris disks (Thebault & Augereau 
2007 ; Krivov et al. 2006 ), can be thoroughly investigated. In ad- 
dition, ALMA will enable detection and imaging of many faint 
disks, and permit detailed SEP modelling of the very coldest 
Kuiper-Belt systems. 

5.3. Dust Masses and Characteristic Radial Distances 

For an optically thin disk at submm wavelengths we can make 
an estimate of the mass of large and cool dust grains (which con- 
tain most of the dust mass), providing a lower limit on the total 
dust mass required to reproduce the measured flux ( Hildebrand 
[1983) : 

F v d 2 
M dust = — - — — -, (2) 

K v hl v (l dust) 

where d is the distance to the source, and the integrated flux den- 
sity F v at 870 pm is assumed to lie on the Rayleigh-Jeans tail of 
the SEP, thus B v appropriately expressed as B v = 2v 2 kT/c 2 . 
The opacity (or dust mass absorption coefficient), k y at submm 



wavelengths is currently unknown due to uncertainties in the 
composition and structure of dust grains, but based on results 
from Praine & Lee ( 1984 1 the 870/vm opacity for micron-sized 
spherical silicate grains should be ~0.5cm 2 g _1 , and approxi- 
mately four times as high for graphite grains at a temperature of 
100 K. Allowing for porous grains would significantly increase 
the opacity (see Stognienko et al. 1995, and references therein), 
but we ado pt K v —2cm £ g~ l (as in Liseau et al.||2003[ |Nilsson| 
et al.||2009), similar to the commonly chosen 1.7 cm g -1 (e.g. 
Zuckerman & Becklin|1993||Pent et al.|2000||Naj"ita & Williams 



2005 ). The distance d to the source and the derived dust temper- 
ature r^ust are two additional parameters with large errors that 
amplify the inaccuracy of the mass determination. Calculated 
Md us i is given in TablePJ] with errors based only on the uncer- 
tainty in the integrated flux density. 



In general, the estimated dust masses fall in the mass range of 
a few moon masses which is typical for submm detected debris 
disks at current sensitivities. One exception is HP 95086 with a 
dust mass in excess of 76Mmooii assuming a second disk com- 
ponent at the upper 45 K temperature limit. 



Another disk parameter that can be obtained from the best-fit 
^dust andyS is the characteristic radial dust distance Rd ust , which 
can be found by assuming thermal equilibrium of dust grains 
(see e.g. Emerson 1988), i.e. the emitted (modified blackbody) 
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In <|3j we made the variable substitution x = (hv)/(kTd ust ). The 
integral can be analytically expressed yielding a characteristic 
radial dust distance given by 



"dust — 



L t {l-A)c 2 ^ 



1/2 



32n 2 h (kT du Jh) 4+ ? £, (4 + ft) T (4 + B) ) 



(5) 



where L» is the luminosity of the star, and A is the albedo of the 
dust grains. g(s) = Y.™ =] rT s is the Riemann zeta function and 
T(n) is the Gamma function. Here we also have to assume a vo 
and corresponding grain size a for which the adopted power-law 
dependence of the opacity is valid. According to the results from 
|Draine & Lee| ( [T984| the power-law relation holds for Aq<,3/mi, 
and we choose a = 1.0 /xm in our calculation. The albedos of 
micron-sized dust grains are generally below 0.1, as shown by 
modelling ( |Shen et al.|2009) and observatio ns of cometary and 



zodiacal dust (Lasue & Levasseur-Regourd 2006 Reach et al 
|2003| l, thus we use A = 0, and note that an albedo of almost 
0.20 would be required to reduce R dust with 10%. The results are 
presented in the last column of Table|2] With the exception of 
HD 98800B, for which we derive a characteristic dust distance 
of ~2AU in accordance with e.g. |Furlan et al.| ( |2007) ; |Akeson| 
|et al.| ( |2007| l; |Prato et aL] ( |2001|i and consistent with the smallest 
stable circumbinary orbit (Holman & Wiegert 1999), all detected 
disks have R dust comparable to Kuiper-Belt distances or beyond. 
The inferred sizes of the detected disks is another reason for 
our hesitation in regarding the extended flux density distribu- 
tion around HD 109573 and HD 141569 as resolved disks. At 
distances of 67 pc and 99 pc, and sizes of 77 AU and 110AU, 
respectively, their angular size in comparison with some other 
disks (e.g. the 170 AU disk around HD 170773, 36 AU away) 
would be smaller (assuming similar grain properties and size dis- 
tributions), suggesting these other disks should also have been 
resolved in our observations. 



5.4. Fractional Dust Luminosities 

Since stellar distances, dust opacities and temperature distribu- 
tions are uncertain, the error in the determined Md ust is actu- 
ally larger than the ones given in Table|2] (which was derived 
from measured flux-density errors), but difficult to quantify. A 
more independent estimate of the amount of dust in the sys- 
tem can be made by calculating the fractional dust luminosity, 
/dust = ^dust/it, from the fitted SEDs. In Table|2]we have listed 
fractional dust luminosities found from our SED modelling of 
detected sources, and values from the literature for undetected 
sources. Again HD 98800 stands out, with an exceptionally large 
/dusb attributed to hot dust belts on the order of ~ 1 AU from the 
binary HD 98800B ( |Furlan et al.|2007) . 

Since all observed stars are nearby the effects of interstellar 
extinction should be be negligible, which is also confirmed by 
B - V colour excesses below ~0. 1 . The only exception is again 
HD 152404 (which is however not modelled) with E B -v - 0.28 . 



This could be both due to its distance (145 pc) and its T Tauri na- 
ture. The general consequence of uncorrected reddening would 
be an overestimation of the fractional dust luminosity, and has to 
be checked for in sources outside the Local Bubble. 



5.5. Temporal Evolution of Debris Disks 

Both theoretical and previous observational work have indicated 
that the amount of dust in debris disks, on the average, will de- 
crease from ~10Myr as the disk ages. This is to be expected, 
due to the combined effect of collisional grinding producing 
small grains that are blown out by radiation pressure, Poynting- 
Robertson drag, sublimation and photosputtering of ices (e.g. 
Wyatt et al.|2007| |Wyatt|2008| l, an d is the general behavior c on 



with a 



firmed in debri s disk surveys (e.gJ N ajita & Will iams 2005 Su 
|et al.|2006[ >. |Spangler et~al| ( |2001) found / dust ocf " 
from ISO observations of clusters, but these results were not con 
firmed by later studies that found significantly lower a, e.g. 



Liu 



leTaTlpXHl i derived a = 0.5-1.0 and |Su et al.| ( [2006l > a = 0.6, 
from Spitzer and JCMT observations, respectively. The main un- 
certainty lies in the inaccuracy of current age estimates of indi- 
vidual stars, as discussed in Section|2j and possible variations 



Nilsson 



in disk evolution depending on stellar spectral types 
|et al. ( 2009[ > compared the average fractional dust luminosity 
observed in the ~12Myr B Pictoris Moving Group with the 
~100Myr Pleiades cluster (Greaves et al. 2009) and obtained 
0.8. Interestingly, these values are higher than those pre- 



a > 



dieted from models of collisional dust evolution, which predict 
a slower decline, with a as low as 0.4 ( |Wyatt et al.|200 7 ; Lohne 
|et al.|2008] >. In Fig.|3]we have plotted /d ust versus age for the ob- 
served sample of stars in this study, with the diameter of the cir- 
cular marker scaled linearly with stellar spectral type (increasing 
for earlier spectral types). Due to the large age uncertainties (rep- 
resented by horizontal error bars) we refrain from making a fit, 
and instead only display dashed lines showing the expected de- 
cline for three different values of a. We conclude that the tempo- 
ral evolution of the fractional dust luminosity determined from 
our observations is in agreement with Nilsson et al. (2009 ), and 
/dust could perhaps decline as fast as the t~ L found by Sp angler| 
et al. (2001), although a < 0.8 can not be excluded due to large 
error bars and small stellar sample. This trend will be better con- 
strained after the completion of the ongoing APEX/LABOCA 
Large Programme. 

Note that this issue is an important one because the value 
of a should be an indication of the dominant process for dust 
removal: a > 1 values are in principle only possible for systems 
at steady-state where Poynting-Robertson drag dominates, while 
a < 1 corresponds to collision dominated systems at steady state 
(Dominik & Decin 2003 ). However, the problem could be more 
complex if the collisional evolution of the system is controlled 
by self-stirring, triggered by the formation of large, 1000- 
2000km sized embryos. In this case, the system can brighten 
from an early low state before later reaching the steady state col- 
lisional decline ( |Kenyon & Bromley|2008l [20041 |Wyatt|2008) . 
An interesting characteristic of self-stirred disks is that the peak 
of dust production, which coincides with the formation of plane- 
tary embryos, should be reached at different times depending on 
the radial distance to the star. This dust production peak should 
thus propagate outward with time, leading to a progressive in- 
crease of the observed characteristic dust distance R dust with age 
as f 1 / 3 . We have tried to identify this evolutionary trend by plot- 
ting /?d us , of our detected submm disks versus age (in Fig. HI. 
In contrast to |Najita & Williams ( 2005 1, who found no appar- 
ent correlation for a sample of 14 submm detected disks, we do 
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Fig. 3. Fractional dust luminosity plotted against stellar age for 
the observed sample of stars. For detected sources our derived 
values have been used, while values from the literature (where 
available) have been used for undetected ones (see references 
in TableEl. Age references are given in Table[TJ The size of 
the circular markers represent the stellar spectral type (diameter 
increase linearly with earlier spectral type). Dashed lines have 
been inserted for comparison with the temporal evolution found 
from previous observations and models. 



identify an increase in characteristic radial dust distance with 
increasing age, although it is again difficult to make any firm 
conclusion regarding the power-law exponent. These tentative 
results will also have to be confirmed with additional detections 
in ongoing observations. 




10' 

Age (Myr) 



Fig. 4. Characteristic radial dust distance for detected submm 
disks plotted against stellar age. Age references can be found in 
TablefT] The dashed line shows the expected slope for the tempo- 
ral evolution of dust in a planetesimal disks dominated by self- 
stirring ( |Kenyon & Bromley|2008"l [2004] ) . 



For both the fractional dust luminosity and the characteristic 
radial dust distance there seems to be no obvious dependence on 
stellar spectral type. 

6. Conclusions 

The most important results from our 870 //m observations of 22 
exo-Kuiper-Belt candidates can be summarised as follows: 

• Out of the observed sample of far-IR excess stars we de- 
tected ~45% (10 out of 22) with at least a 3<x significance. 
Five of the detected debris disks (HD 95086, HD 131835, 
HD 161868, HD 170773, and HD 207129) have previously 
never been seen at submm wavelengths. These findings in- 
crease the currently known exo-Kuiper-Belts (inferred from 
cold extended dust disks detected in the submm) from 29 
to 34. The sensitivity of current submm observations with 
LABOCA at APEX and SCUBA at JCMT are reaching the 
background extragalactic confusion limit, as evident by the 
abundance of significant off-center flux density peaks in ac- 
quired submm maps. Future ALMA observations will be 
essential to discriminate between resolved disks and back- 
ground submm galaxies. 

• We perform SED modelling of available photometric data 
using a simple modified blackbody to make a x 2 minimiza- 
tion fit to the disk excess emission in order to determine 
the dust temperature rd ust and power-law exponent /? of the 
opacity law. Resulting /3 parameters are all between 0. 1 and 
0.8, suggesting grains significantly larger than those in the 
ISM, in agreement with previous results for debris disks (e.g. 
Najita & Williams 2005 ; Nilsson et al |2009| |Roccatagliata| 



et al. 



2009). However, it should be noted that the presence 



of co der secondary dust components in addition to the fitted 
disk could be responsible for some abnormally low [5 values. 
The measured 870 fim flux of HD 95086 undeniably requires 
a second disk component. 

Besides finding the dust mass of detected sources (together 
with upper limits on undetected ones) derived from the in- 
tegrated 870 //m flux density, we use the best-fit SED pa- 
rameters to calculate the fractional dust luminosity /d ust and 
the characteristic radial dust distance 7?dust (for 1/mi sized 
grains) of the disks. We present plots of the temporal evo- 
lution of /dust and 7?dust based on these values together with 
previously estimated stellar ages, and find that our results are 
consistent with a /d ust °c t~ a , with a ~ 0.8-2.0, as well as ten- 
tatively with the Rdust <* t l ^ relation predicted in debris disks 
collisional evolution models by Kenyon & Bromley (2004). 
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Appendix A: Spatially Resolved Disks or Chance 
Alignment of Background Galaxies? 

Two sources appear spatially resolved in our 870 pm maps, both 
of them are multiple systems which have been previously im- 
aged at shorter wavelengths. Here we compare the resolved 
structures at submm wavelengths with that observed in scattered 
optical light and IR emission, and discuss the indications for and 
against them actually being resolved disks. The position angle 
(PA) was found by fitting a 2-D Gaussian profile to the extended 
submm emission. 



A. 1. HD 1 09573 (HR 4796) 

The highly inclined debris disk around HD 109573 A (also 



known as HR 4796A) was first imaged in the mid-IR by|Koerner 



etal.| ( |1998) and |Jayawardhana eta l. ( 1998), followed by optical 
and near-IR coronagraphic imaging by Sc hneider et al.| ([l999 
2009). These observations revealed a 17 AU wide dust ring with 
~70AU radius surrounding the primary A0V star. In our map 
(Fig. |l(f)[ ), the peak flux density of 19.1 ±3.1 mJy/beam is cen- 
tered on the position of the primary, with the secondary M2.5 
star located 7/7 (~2 pixels) to the southwest. Surprisingly, the 
PA is 129°, compared to 27° for the ring imaged by Schneider 
et al.| ( |2009[ l, i.e. a nearly orthogonally oriented extension. 
Although the accuracy of our determined PA is hard to estimate 
due to the influence of possible lower-sigma noise features, the 
flux distribution is undoubtedly extended in the southeast direc- 
tion, maybe with a warp toward the East. The 3cr contour reaches 
>2000AU projected radial distance on that side of the central 
star. However, a chance alignment with a bright background 
galaxy to the southeast is probable, for several reasons. Firstly, 
an outer disk orthogonally oriented to the inner one seems highly 
implausible or even unphysical. In addition, a data reduction em- 
ploying only half-beam smoothing seem to imply a separate 3<x 
peak in the southeast flux extension. A second galaxy to explain 



the small northwest extension is conceivable, but it could also be 
a noise feature. The formal likelihood of a source falling within 
some 25" distance from the star and making it appear extended 

is about 4%. 

A.2. HD 141569 

This triple system has been extensively studied since the Hubble 
Space Telescope (HST) discovery of a large disk surrounding 



the primary AOVe star HD 14 1569 A (Weinber ger et al.|[1999 
Augereau et al.|[T999a|l. The disk h as been imaged in thermal 



mid-IR emission (Fisher et a l.|2000[ l, but up until now never re- 
solved at longer wavelengths. We find a clear asymmetric elon- 
gation of the disk with PA of 10° to the southwest of the source- 
centered 10.3±2.3mJy/beam peak flux (Fig. l(h)| i. This can be 
compared to the two most detailed studies of the disk structure so 
far, employing coronagraphic HST ACS and STIS observations, 
respectively (Clampin et al. 2003 Mouillet et al. 2001)1, which 



revealed an inner (<175 AU) clearing surrounded by two ring- 
like structures - a thin (~50AU) belt at ~200AU and a wider 
(~100AU) belt at ~350AU - oriented with a semi-major axis 
roughly in the North-South direction and connected by faint spi- 
ral structures. Clampin et al. (2003 1 also found additional spi- 
ral arcs extending from the northeast (out to ~1200 AU) and the 
southwest (toward the binary HD 141569BC, located ~8" north- 
west of HD 141569A. The reason for the highly asymmetric 
brightness distribution in the optical, with respect both to semi- 
major and semi-minor axis, could be a non-axisymmetric distri- 
bution of grains and/or anisotropic scattering by grains (Mouillet 
et al. 2001). Suspected gravitational perturbation by a massive 
body in the disk, and effects from the binary companion has 
been extensively stud i ed by e .g. Augereau & Papaloizou (2004); 
Quillen et al] ([2005); |Wyatt| ( |2005[ ), while |Ardila et"aL]p005) 



and |Reche et al.| ( |2009| l thoroughly investigated the case of a 
gravitationally unbound binary passing by (with and without 
planets in the disk), but did not manage to explain all observed 
features. It is possible that we are seeing a colder dust popula- 
tion on very eccentric orbits, perhaps shaped by the dynamical 
interaction with the binary (either in a flyby or a triple system 
scenario), however, dynamical modelling of such a scenario is 
beyond the scope of the present paper. The extreme brightness 
asymmetry with respect to the semi-minor axis, on the other 
hand, does not make us confident that the extended emission 
originates from the dust disk, and could instead be attributed 
to a background galaxy located southwest of HD 141569. Also 
in this case a re-reduction applying only half-beam smoothing 
seem to hint at a separate 3cr peak, only 8" to the southeast of 
the star (which would be a 0.5% chance alignment). 

Appendix B: List of submm detected 
exo-Kuiper-Belts 
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Table B.l. Submm and mm photometry of all debris disks that have been detected at submm and longer wavelengths. 



Star 


Other Name 


Gum) 


(Jy) 


Reference 


HD105 




870 


10.7 ± 5.9° 


1 


HD377 




1200 


4.0 ± 1.0 


11 


HD 8907 




1200 


3.2 ±0.9 


11 


HD 14055 


y Tri 


850 


5.5 ± 1.8 


3 


HD 15115 




850 


4.9 ± 1.6 


3 


HD 17206 


t 1 Eri 


1300 


20.7 ±3.9 


6 


HD 21997 




870 


8.3 ±2.3 


3 






870 


17.6 ±8.0 


1 


HD 22049 


eEri 


450 


225 ± 10 


12 






850 


40.0 ±1.5 


12 






1300 


24.2 ±3.4 


6 


HD 30447 




870 


6.9 ±5.0" 


1 


HD 32297 




1300 


5.1 ± 1.1 


10 


HD 38393 




850 


2.4 ± 1.0 


12 


HD 39060 


ySPic 


850 


58.3 ±6.5 


4 






870 


63.6 ±6.7 


2 






1200 


24.3 ± 3.0 


5 






1300 


24.9 ± 2.6 


6 


HD 39944 




1350 


3.7 ±0.9 


14 


HD 48682 


56 Aur 


850 


5.5 ± 1.1 


12 


HD 61005 




350 


95 ± 12 


11 


HD 95086 




870 


41.3 ± 18.4 


1 


HD 98800* 


TVCrt 


870 


33.6 ±8.4 


1 






1350 


36.8 ±4.2 


14 






2000 


24.8 ± 3.4 


14 


HD 104860 




350 


50.1 ±9.3 


11 






1200 


4.4± 1.1 


11 


HD 107146 




350 


319±6 


11 






450 


130 ±40 


7 






850 


20 ±4 


7 


HD 109085 


77 Cor 


850 


7.5 ± 1.2 


12 


HD 109573 


HR4796 


450 


180 ± 150 


12 






850 


19.1 ±3.4 


8, 12 






870 


21.5 ±6.6 


1 


HD 123160 




850 


13 ±4.3 


14 






1350 


4.7 ±0.9 


14 


HD 128167 


cr Boo 


850 


6.2 ± 1.7 


12 


HD 131835 




870 


8.5 ±4.4 


1 


HD 141569* 




450 


66 ±23 


12 






850 


14 ±1 


12 






870 


12.6 ±4.6 


1 






1300 


9.39 ± 3.72 


13 






1350 


5.3± 1.1 


14 


HD 152404* 


AKSco 


450 


242 ± 32 


15 






800 


83 ±9 


15 






870 


42.9 ±9.8 


1 






1100 


36 ±4 


15 


HD 161868 


7 Oph 


870 


12.8 ±5.2 


1 


HD 170773 




870 


18.0 ±5.4 


1 


HD 172167 


Vega 


850 


45.7 ±5.4 


4, 12 






1300 


11.4 ±1.7 


16 


HD 181327 




870 


51.7 ±6.2 


2 


HD 191089 




350 


54 ± 15 


11 


HD 195627 




870 


13.0 ±7.1" 


1 


HD 197481 


AUMic 


850 


14.4 ±1.8 


9 


HD 207129 




870 


5.1 ±2.7 


1 


HD 216956 


Fomalhaut 


450 


595 ± 35 


12 






850 


97 ±5 


12 






850 


81 ±7.2 


4 


HD 218396 


V342 Peg 


850 


10.3 ± 1.8 


3 


HIP 23200 


VI 005 Ori 


850 


4.8 ± 1.2 


9 



References. (1) This paper; (2) Nilsson et al. (2009); (3) Williams & Andrews 
|eTaL l(fT99i); (7)[Williams et al.|(|2004[ 



2006); (4) Holland et al. ( 1998); (5) |Liseau et af](i 2003 >; (6)jCEm 
2004); (10) Maness et al. (2008); (11) Roccatagliata et al.|([2009" 



(8)|Greaves et al.H2000>; (9)|Liu et al.| , 
(12) |Sheret et al.] ( |2004| ); (13) |Walker & Butner|11995] >r(14) |Sylvester et al.| l |200T] >rtl5) |Ten"sen et al.| < |199^ >ra6) |WUn"er et al.H2002'] l 

Notes. (o) Marginally detected. {b) Debris disk status disputed. 
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